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ABSTRACT 

This  report  discusses  the  design  and  impletaentation  of  FORTRAN  subrou> 
tines  to  add  the  capabilities  of  nonseparability  testing  and  pivotal  decomposi* 
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1.  Introduction 

)  PolyChain  is  a  portable  FORTRAN  program  for  evaluating  the  reliability  of  a  K-tenninal 
netvk'ork  via  polygon-to-chain  reductions  and  the  factoring  algorithm.  The  first  version  of 
PolyChain  allows  the  evaluation  of  the  K-terminal  reliability  for  series-parallel  graphs  [1]:^ 
The  algorithm  implemented  is  a  linear  time  algorithm  introduced  in  1982  by  Satyanarayana 
and  Wood  12, 3]^ 

This  report  discusses  the  design  and  implementation  of  two  features  recommended  in^|]|^that 
enable  PolyChain, to  treat  a  larger  class  of  problems.  The  algorithm  of  Satyanarayana  and 
Wood  has  a  constraint  on  the  topology  of  the  input  network,  requiring  it  to  be  nonseparable. 
The  original  version  of  PolyChain  does  not  test  for  this  requirement.  One  of  the  features 
added  to  PolyChain  discussed  in  this  report  is  the  implementation  of  a  routine  to  check  for 
this  condition. 

The  algorithm  of  Satyanarayana  and  Wood  computes  the  reliability  of  a  network  with  an 
underlying  series-parallel  structure.  When  the  input  netwoik  is  not  totally  reducible  an  exten¬ 
sion  to  the  algorithm  is  required  to  obuin  the  reduced  network.  The  second  feature  discussed 


in  this  report  is  the  implementation  of  a  factoring  algorithm  incorporated  to  PolyChain  to 
insure  the  evaluation  of  the  K-terminal  network  reliability  for  both  series-parallel  reducible 
and  irreducible  networks.  — — 

Section  2  briefly  presents  some  theoretical  results  of  polygon-to-chain  reductions.  In  section 
3,  the  algorithm  implemented  for  nonseparability  testing  is  presented.  Section  4  briefly 
discusses  the  factoring  algorithm.  A  system  manual  is  presented  in  section  S  describing  the 
implementation  of  both  algorithms  in  FORTRAN,  and  a  user  manual  is  presented  in  section 
6.  The  code’s  performance  is  illustrated  in  section  7  through  the  testing  of  several  networks. 
Conclusions  and  recommendations  are  made  in  section  8. 

2.  Series-Parallel  Graphs  and  Polygon-to-Chain  Reductions 

In  this  section  a  brief  discussion  of  series-parallel  graphs  and  polygon-to-chain  reductions  is 
presented.  For  a  complete  discussion  see  [2,3]. 

Throughout  this  report  we  consider  an  undirected  graph  G={VJE),  where  V  is  the  set  of  ver¬ 
tices  and  E  the  set  of  edges  of  G.  A  connected  graph  G=(K,£)  is  said  to  be  separable  if 
there  exists  a  vertex  v,  called  the  separation  vertex,  such  that  its  removal  from  the  graph 
disconnects  the  graph.  When  a  graph  has  no  separation  vertex,  it  is  called  nonseparable .  The 
induced  nonseparable  subgraphs  of  a  separable  graph  G  are  called  nonseparable  components 
ofG. 

Let  G=(V^)  be  a  nonseparable  graph.  Vertices  are  assumed  to  be  perfectly  reliable,  and 
edges  may  fail,  independently  of  each  other,  with  known  probabilities.  The  edge  reliability 
for  edge  e,  is  p,  ,  and  the  edge-failure  probability  is  9,  «I-Pi.  htx  K  ^V,  |Ar  |  ^  2  be  a 
specified  set  of  vertices.  Vertices  in  K  will  be  referred  to  as  K -vertices.  Gg  is  graph  G  with 
K  specified.  The  K-terminal  reliability  of  Gk,  R{Gk),  is  the  probability  that  all  K- vertices 
in  Gic  are  connected  by  working  edges. 


The  size  of  graph  Ca,  i  e.  |  ^'(Ga)  I  +  I  ^^(Ga)  I .  can  be  reduced  by  applying 
reliability  -preserving  reductions.  The  application  of  reliability-preserving  reductions  to  Ga 
renders  a  graph  G'a  such  that  /?(Ga  )  =  fi/?(G  a  ),  where  Q  is  a  multiplicative  factor  that 
depends  on  the  reductions  applied. 

Three  types  of  reliability-preserving  reductions  will  be  referred  to  as  simple  reductions:  paral¬ 
lel  reduction,  series  reduction,  and  degree-2  reduction.  In  parallel  reduction,  parallel  edges 
ea=(.x,y)  and  eb={x,y)  are  replaced  by  a  single  edge  ec’‘{x,y)  with  edge  probability 
Pc  =  ^-QaQh-  Jo  series  reduction,  edges  ea={x,y)  and  ei,=(y,z)  are  replaced  by  a  single  edge 
ec=(x,z)  with  edge  probability  Pc=PaPb-  both  series  and  parallel  reductions  the  multipli¬ 
cative  factor  n  has  value  1,  and  K'=K.  In  degree-2  reduction  edges  ea^{x,y)  and 
x,y,z€K,  are  replaced  by  edge  ec={x,z)  with  Pc=PaPb!(^-da<lb\  and 

K'=K  -  y. 


Replacing  a  pair  of  series  (parallel)  edges  by  a  single  edge  is  called  a  series  (parallel)  replace¬ 
ment.  A  replacement,  as  opposed  to  a  reduction,  does  not  involve  probabilities  or  a  set  of 
distinguished  nodes  associated  with  the  graph. 


A  nonseparable  series  -parallel  graph  is  a  graph  that  can  be  reduced  to  a  single  edge  by  suc¬ 
cessive  series  and  parallel  replacements.  If  the  graph  is  separable,  it  is  series-parallel  if  it  can 
be  reduced  to  a  tree  after  all  possible  series  and  parallel  replacements  are  performed.  A  non¬ 
separable  series-parallel  graph  Ga  is  termed  s-p  reducible  if  it  can  be  reduced  to  a  single 
edge  by  successive  simple  reductions.  A  graph  Ga  is  s-p  irreducible  if  it  is  not  s-p  reduci¬ 
ble. 


A  chain  is  an  alternating  sequence  of  distina  vertices  and  edges,  such  that  the  internal  ver¬ 
tices  are  all  of  degree  2  and  end  vertices  are  of  degree  greater  than  2.  A  chain  must  contain  at 
least  one  edge  and  two  end  vertices.  A  polygon  is  a  cycle  such  that  exactly  two  vertices  of  the 
cycle  are  of  degree  greater  than  2. 
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A  set  of  reliability-preserving  reductions  introduced  by  Satyanarayana  and  Wood  [2,3], 
replaces  a  polygon  with  a  chain.  These  reductions  are  called  polygon -to -chain  reductions. 
It  is  shown  in  [2,3],  that  every  series-parallel  graph  is  reducible,  irrespective  of  the  vertices 
chosen  to  be  in  K,  with  the  use  of  simple  reductions  and  polygon-to-chain  reductions.  Mak¬ 
ing  use  of  these  two  types  of  reliability-preserving  reductions,  a  linear  time  algorithm  to 
evaluate  R(Gtc)  for  a  series-parallel  graph  with  any  chosen  set  K  is  presented  in  [2,3]. 

PolyChain  is  a  direct  implementation  of  that  algorithm  utilizing  an  extension  so  that  a 
reduced  network  can  be  obtained  when  the  graph  is  s-p  irreducible.  When  a  reduced  graph  is 
generated  the  factoring  algorithm  is  applied  to  find  the  reliability  of  the  reduced  network. 
This  way.  Polychain  can  evaluate  the  K-terminal  reliability  of  general  nonseparable  networks. 

3.  Nonseparability  Testing 

A  depth-first-search  based  algorithm,  having  time  complexity  CHIEI),  exists  to  detect  separat¬ 
ing  vertices.  This  algorithm  is  implemented  in  PolyChain  and  is  presented  below  as 
described  in  [9]. 


Assume  that  |  |  >  1,  and  .r  is  the  vertex  in  which  we  start  the  search. 

(1)  Mark  all  edges  'unused'.  Empty  the  stack  5. 

For  every  v  €  F  let  k(v)  •*-  0.  Let  i  0  and  v  *-  s. 

(2)  i  -^1  +  1,  k(v)  i,  L(v)  /  and  put  v  on  S. 


All  the  vertices  on  5  down  to  and  including  v  are  now  removed 
from  S;  this  set,  with  /(v),  forms  a  nonseparable  component. 

(8)  V  /(v)  and  go  to  Step  (3). 

(9)  All  vertices  on  5  down  to  and  including  v  are  now  removed 
from  S ;  tncy  form  with  s  a  nonseparable  component . 

(10)  If  5  has  no  unused  incident  edges  then  halt. 

(11)  Vertex  5  is  a  separating  vertex.  Let  v  ^  s  and  go  to  Step  (4). 

4.  Factoring  Algorithm 

As  already  mentioned,  when  the  input  network  has  no  underlying  series-parallel  structure,  the 
polygon-to-chain  algorithm  generates  a  reduced  network  but  does  not  compute  the  network's 
K-terminal  reliability.  We  will  discuss  the  implementation  of  the  factoring  algorithm  incor¬ 
porated  to  PolyChain  to  calculate  the  reliability  of  the  reduced  network  generated  when  the 
input  network  is  s-p  irreducible. 

4.1.  The  Algorithm 

The  K-terminal  reliability,  R^iG),  of  a  graph  G  can  be  computed  by  repeated  applications  of 
the  following  decomposition, 

Rk{G)=p,R(G,)H\-p,)R{G.,) 

where  Gf  is  the  graph  obtained  from  G  by  considering  that  edge  e  is  working  and  G.,  is  the 
graph  obtained  from  G  when  edge  e  is  not  working.  Hence,  G,  and  G.,  are  obtained  by 
respectively  contracting  and  deleting  edge  e  in  G. 

After  each  application  of  this  decomposition,  simple  reductions  are  performed.  If  the  gen¬ 
erated  subgraph  is  not  totally  reduced  a  new  edge  is  then  selected  and  the  decomposition 


reapplied. 


The  following  schente  describes,  in  recursive  form,  the  factoring  algorithm. 

factor  (G) 
reduce (G) 
select  edge  e  to  pivot 
factor  (Gr) 
factor  {G-e) 
end 

The  use  of  the  factoring  algorithm  generates  a  binary  computational  tree  whose  root  node  is 
the  original  graph  and  each  other  node  is  a  sul^raph.  Without  the  application  of  simple 
reductions  after  each  edge  selection,  the  binary  structure  would  contain  2 1  ^ '  leaves,  which  is 
equivalent  to  the  enumeration  of  all  possible  states  of  G.  Notice  that  the  scheme  given  above 
traverses  the  binary  computational  tree  using  a  preorder  enumeration.  An  example  of  a 
preorder  traversal  applied  to  a  tree  is  given  in  figure  1. 

For  a  complete  discussion  of  the  factoring  algorithm  see  [4,5]. 

PREORDER:  A  B  D  E  G  C  F 

Figure  1 


4.2.  Edge  Selection 

In  this  section,  results  on  optimal  edge  selection  of  Satyanarayana  and  Giang  [10]  are 
reviewed.  Satyanarayana  and  Chang  show  that  there  exists  an  edge  selection  that  yields  the 
optimal  binary  structure,  that  is,  a  binary  tree  with  the  minimal  number  of  leaves.  They  call 
such  an  edge  selection  the  optimal  edge  selection.  They  also  show  that  the  number  of  leaves 


of  the  optimal  binary  tree  is  equal  to  the  domination,  Dic(G),  of  G.  A  K-tree  is  a  tree  of  G 
covering  all  K-vertices  and  having  its  pendant  vertices  in  An  irrelevant  edge  is  one  that 
lies  on  no  K-tree.  An  edge  selection  is  optimal  if  and  only  if  every  reduced  graph  generated 
has  no  irrelevant  edge.  Hence,  a  fast  edge  selection  strategy  that  avoids  creating  subgraphs 
with  irrelevant  edges  is  desired. 

A  graph  G,  with  respect  to  some  set  K,  is  termed  a  K -graph  if  every  edge  of  G  is  relevant, 
i.e.,  if  every  edge  of  G  is  in  some  K-tree  of  G.  Satyanarayana  and  Chang  prove  that  for  a  K- 
terminal  irreducible  graph  G  with  domination  Z)jt^(G)>l,  there  exists  an  edge  such  that  Ge 
and  G.f  are  both  K-graphs.  They  further  show,  that  an  edge  satisfying  the  property  men¬ 
tioned  above,  can  be  found  in  0(  |  £  |  +  |  K  | )  operations  using  techniques  based  on  depth- 
first-search  [11,12]. 

In  the  version  of  the  factoring  algorithm  implemented  in  PolyChain  the  optimal  edge  selec¬ 
tion  strategy  is  used.  At  each  iteration  an  edge  whose  removal  does  not  disconnect  the  graph 
is  chosen.  Then,  each  of  the  subgraphs  generated  by  pivoting  on  the  selected  edge  is  checked 
for  irrelevant  edges.  Notice  that  a  graph  is  a  K-grapb  if  and  only  if  each  one  of  its  pendant 
nonseparable  components  has  at  least  one  distinguished  node.  Therefore,  the  algorithm 
checks  if  the  subgraphs  are  nonseparable.  In  the  case  a  graph  is  separable,  the  algorithm 
checks  if  all  of  its  pendant  nonseparable  components  have  at  least  one  distinguished  node  in 
it.  If  a  pendant  nonseparable  component  not  having  any  distinguished  node  exists  in  either 
one  of  the  subgraphs  generated  by  the  edge  selected,  the  current  edge  is  discarded  and  another 
edge  is  selected  to  replace  it.  Then,  the  checking  procedure  starts  all  over  again  considering, 
now,  the  new  edge  selected.  Since  a  graph  having  only  one  distinguished  node  is  not  a  K- 
graph  the  algorithm  avoids  creating  such  graphs.  This  edge  selection  procedure  is  0(  |  £  |  ^). 


■  v^L  i v^y  v^if^v* miv* < ■  i.  * *y,»i 


1 .  • 


-8- 


5.  System  Manual 


>.-»v 


Liiii 


In  this  section  the  code  is  briefly  described. 


i- 


: 


► 


5.1.  Programming 

All  subroutines  in  this  version  on  PolyChain  are  written  in  Fortran  77.  Only  I/O  related  code 
is  system  dependent.  The  input  network  for  the  factoring  algorithm  is  the  output  of  the  origi¬ 
nal  version  of  PolyChain  when  its  input  network  is  not  totally  reduced. 

5.2.  Data  Structures 

PolyChain  uses  an  efficient  network  representation  using  linked  list  data  structures,  [6,7,8]. 
Each  vertex  has  a  list  of  adjacent  vertices,  which  not  only  indicates  which  vertices  are  adja¬ 
cent  to  it,  but  also  provides  information  whether  the  vertex  belongs  to  set  K.  For  every  ele¬ 
ment  of  the  list,  there  is  a  pointer  giving  the  address  of  the  information  about  the  edge.  Fig¬ 
ure  3  illustrates  this  multilist  structure  for  the  network  given  in  6gure  2. 

The  routines  incorporated  to  PolyChain  use  a  few  other  data  structures  in  addition  to  the 
ones  used  by  PolyChain .  For  nonseparability  testing  a  stack  data  structure  is  used.  For  the 
factoring  algorithm  another  stack  is  used  to  provide  information  about  the  computational 
binary  tree. 


♦*.  • 


%"• j'" V-  **  •-**  '*’■  j*^-i** '  •' •  ■ 


Figure  3 


5.3.  Data  Structure  Implementation 

Next,  we  describe  the  FORTRAN  arrays  used  to  implement  the  data  structure  of  the  two  rou> 
tines. 

5.3.1.  Nonseparability  Test 

As  already  mentioned  a  stack  is  used  for  producing  the  vertices  of  the  component.  The  ver¬ 
tices  are  stored  in  the  stack  in  the  order  that  they  are  discovered.  When  a  separating  vertex  u 
is  discovered,  we  read  off  all  the  vertices  from  the  top  of  the  stack  down  to  a  node  specified 
by  the  algorithm.  All  these  edges  plus  the  separating  vertex  u  constitute  the  component. 

53.2.  Factoring  Algorithm 

Five  arrays  are  used  to  implement  the  stack  that  stores  information  about  the  computational 
binary  tree.  EDGEV1(*)  and  EDGEV2(*)  contain  the  stack  of  vertices  corresponding  to  the 
selected  edge.  D1RECT(*)  contains  information  about  the  branching.  If  DIR£CT(*)  is  1,  the 
selected  edge  is  working.  If  it  is  -1,  the  selected  edge  is  not  working.  EPROB(*)  contains  the 
reliability  of  the  selected  edge,  and  RELB(*)  contains  the  value  of  M  after  all  possible  simple 
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reductions  are  performed.  If  no  degree-2  reduction  is  performed,  the  value  of  M  is  1.  If  a 
degree-2  reduction  is  performed,  the  value  of  M  will  be  updated.  TOP  points  to  the  top  of 
stacks  EDGEV1(»),  EDGEV2(*),  DIRECT(*),  EPROB{*),  and  RELB(*). 

The  implementation  of  the  factoring  algorithm  was  carried  out  in  a  way  to  minimi/.e  core 
usage.  As  already  mentioned,  a  preorder  binary  tree  traversal  algorithm  was  implemented. 
Hence,  after  an  edge  is  selected  we  always  consider  first  the  case  in  which  the  selected  edge  is 
working.  When  the  subgraph  can  be  reduced  by  simple  reductions,  the  algorithm  finds  its 
reliability  and  goes  back  to  its  parent  node  to  continue  branching.  After  leaving  a  node, 
which  is  actually  a  subgraph,  that  was  already  branched  in  both  directions,  the  algorithm 
never  comes  back  to  it.  Hence,  this  subgraph  and  all  subgraphs  beneath  it  do  not  have  to  be 
saved  and  can  therefore  be  deleted.  Since  the  factoring  algorithm  was  not  implemented  in 
recursive  form,  we  have  to  keep  the  information  necessary  for  the  recovery  of  the  subgraph  of 
the  computational  binary  tree. 

To  better  understand  the  method,  suppose  the  computational  binary  tree  is  of  the  form  given 
in  figure  4. 


Figure  4 


The  numbered  nodes  correspond  to  series-parallel  reducible  subgraphs.  The  implemented  fac¬ 
toring  algorithm  finds  the  reliability  of  each  branch  and  then  the  reliability  of  the  original 
graph  as  shown  in  figures  4a,  4b,  4c,  4d,  4e,  4f.  and  4g.  For  example,  if  h,  is  the  reliability  of 
the  branch  leading  to  leaf  t,  the  overall  reliability  of  the  original  graph  (7a  is 
Rk(G)  =  M  m,  where  A/  =  obtained  from  the  polygon  -  to  -chain  reductions.  The 

I  j 

reliability  of  a  leaf,  h,,  is  the  product  of  the  reliabilities  of  the  edges  selected  leading  to  that 
leaf  For  example,  A3  =  Pe\Pe2Pe*Pei- 


^  ^  / 


5.4.  Data  Dictionar}' 


In  this  section,  the  variables  used  in  the  new  subroutines  are  listed,  with  a  brief  description  of 


each  one.  The  variables  marked  with  an  asterisk  are  new  variables,  not  present  in  the  earlier 


version  of  PolyChain .  A  data  dictionary  containing  other  variables  used  in  other  PolyChain 


subroutines  can  be  found  in  [  1  ]. 


ADJVRT(*)  vertex  adjacent  to  vertex  whose  list  it  is  on 


•APTR 


•ATOPFB 


•AUX 


•AUX2 


*AUX3 


•AUXL 


AVSADJ 


CARDE 


CARDV 


•COUNTE 


DATE 


DEG(*) 


•DISTGN 


•DOMINT 


•EDGE 


auxiliary  pointer 


auxiliary  variable 


auxiliary  variable 


auxiliary  variable 


auxiliary  variable 


auxiliary  variable 


pointer  to  beginning  of  list  of  available  space 


•BRIDGE!*)  array  of  edges  whose  removal  disconnect  the  graph 


cardinality  of  set  E 


cardinality  of  set  V 


counter  of  the  number  of  iterations 


degree  of  vertex 


auxiliary  variable 


•DIRECT!*)  direction  of  the  branch  taken 


current  number  of  distinguished  nodes 


domination  of  the  graph 


edge  formed  by  the  given  pair  of  vertices 


•EDGMRK!*)  array  containing  information  about  marked  edges 


•EDGEPV 


edge  that  is  a  candidate  to  be  a  forbidden  edge 


•  i  V- 


•  I 


r  . 


V. 
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•EDGEV1{*)  one  of  the  vertices  of  selected  edge 

•EDGEV2(*)  one  of  the  vertices  of  selected  edge 

EDGPRB(*)  edge  reliability 

•EPROB(*)  array  of  probabilities  of  selected  edges 
*FATHER(*)  preceding  node  in  the  search 

•FTMONSP  indicator  of  whether  or  not  the  output  is  to  be  printed 
•FIRST  starting  node  in  the  search 

•FIRSTV  vertex 

•FORBED(*)  edges  that  are  forbidden  to  be  selected  to  be  pivot 
•FOUND  indicator  of  whether  a  new  edge  exists  to  continue  search 

•FOUND2  indicator  of  whether  a  candidate  degree-2  edge  pair  was  found 

•FOUNDP  indicator  of  whether  parallel  edges  were  found 
•FOUNDS  indicator  of  whether  series  edges  were  found 

•HEAD  auxiliary  pointer 

HOUR  hour 

IN  value  of  FORTRAN  input  file 

lOUT  value  of  FORTRAN  output  file 

•IPTR  pointer 

•K(v)  number  of  vertex  v 

•KGRAPH  indicator  of  whether  or  not  a  graph  is  a  K-graph 

•L(v)  lowpoint  of  v 

•LIMITI  lower  limit 

•LIMITS  upper  limit 


M 

MAXLST 

MCARDE 

MCARDV 

ND2R 

*NUMCMP 

•NUMELM 

•NPIVOT 

•NPR 

NSR 

•OUTPUT 

•POINT 

•PTR 

PTRADJ(^) 

•PTRCMP(^) 

QA 

QB 

•REL 

•RELB(^) 

•SECODV 

•SEPVTX(^) 

•SREL 

•STACK(^) 

•TADJVT(^) 

•TAVSAD 

•TBRIDG(^) 


product  of  all  (1,  see  [3] 

maximum  number  of  elements  in  adjacent  vertices  list 

cardinality  of  set  E  at  start  of  procedure 

cardinality  of  set  V  at  start  of  procedure 

counter  of  degree-2  reductions 

number  of  nonseparable  components 

number  of  elements 

number  of  pivots  performed 

counter  of  parallel  reductions 

counter  of  series  reductions 

network  reliability  after  factoring  algorithm 

pointer 

pointer 

pointer  to  beginning  of  list  of  adjacent  vertices 
pointer  to  beginning  of  list  of  vertices  in  a  component 
failure  probability  of  edge  a 
failure  probability  of  edge  b 
total  reliability 

value  of  M  after  each  simple  reduction 
vertex 

array  of  separable  vertices 
subgraph  reliability 
stack  of  vertices  scanned 
copy  of  current  adjvrt(^) 
copy  of  current  avsadj 
copy  of  current  bridge(^) 


•TCARDE 


copy  of  current  carde 
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•TCARDK 

copy  of  current  cardk 

•TCARDV 

copy  of  current  cardv 

TCPU 

total  epu  time 

•TDEG(*) 

copy  of  current  deg(^) 

•TEDGPB(*) 

copy  of  current  edgprb(^) 

TEST 

key  for  debugging  feature 

•TLNKDW(*) 

copy  of  current  Inkdwn(^) 

*TLNKED(*) 

copy  of  current  lnkedg(^) 

*TLNKUP(*) 

copy  of  current  Inkup(^) 

*TM 

copy  of  current  m 

•TMPCPU 

solution  time 

•TOP 

top  of  chain  stack 

•TOPB 

pointer  to  the  top  of  the  list  of  bridges 

•TOPCMP 

pointer  to  beginning  of  the  component 

•TOPFB 

pointer  to  the  top  of  the  list  of  forbidden  edges 

•TOPS 

top  of  vertex  stack 

•TOPVTX 

pointer  to  the  top  of  the  list  of  separable  vertices 

•TTOPB 

copy  of  current  topb 

•u 

vertex 

•V 

vertex 

•VI 

vertex 

•V2 

vertex 

•VCMP(^) 

vertex 

VCPU 

virtual  epu  time 

•VERTEX 

vertex 

•VERTXl 

vertex 

•VERTX2 


vertex 


VRTX 


vertex 


•XADJVT(*)  adjvrt(*)  at  start  of  factoring  routine 
•XAVSAD  avsadj  at  start  of  factoring  routine 
•XCARDE  carde  at  start  of  factoring  routine 
•XCARDV  cardv  at  start  of  factoring  routine 
•XDEG(*)  deg(*)  at  start  of  factoring  routine 

•XEDGPB(*)  edgprb(*)  at  start  of  factoring  routine 
•XLNKDW(*)  lnkdwn(*)  at  start  of  factoring  routine 
•XLNKED(*)  lnkedg(*)  at  start  of  factoring  routine 

•XLNKUP(*)  lnkup(*)  at  start  of  faaoring  routine 
•XPTADJ(*)  ptradjC*)  at  start  of  factoring  routine 
YEAR  year 

5.5.  COMMON  Blocks 

All  the  COMMON  blocks  used  in  the  new  routines  are  listed  below.  The  COMMON  blocks 
introduced  in  this  new  version  of  PolyChain  are  marked  with  asterisk.  For  a  list  of  all  other 
COMMON  blocks  in  the  code  see  (IJ. 

COMMON/BLKOl/  DEG(MAXVRT) 

COMMON/BLK02/PTRADJ(MAXVRT),ADJVRT(2*MAXEDG), AVSADJ 
COMMON/BLK2 1  /  LNKDWN(2*MAXEDG),LNKUP(2*MAXEDG),LNKEDG(2*MAXEDG) 
COMMON/BLK03/  EDGPRB(MAXEDG),EDGNUM(MAXEDG) 

COMMON/BLK05/  MAXEDG,MAXVRT,MAXLST,MXSTKT,MAXCHN 

COMMON/BLK06/  CARDE,CARDV,CARDK 

COMMON/BLK07/  M 

COMMON/BLK08/  IN.IOUT 

COMMON/BLK3 1  /  MCARDE,MCARDV,MCARDK 


COMMON/BLK32/  DATE.YEAR.HOUR 


•COMMON/BLK40/  K(MAXEDG),L(MAXEDG),STACK(MAXEDG).EDGMRK(MAXEDG),TOPS 
•COMMON/BLK41/PTRCMP(MAXVRT).TOPCMP.NUMCMP,VCMP(MAXEDG) 
•COMMON/BLK50/EDGEV1(MAXEDG),EDGEV2(MAXEDG),D1RECT(MAXEDG) 
•COMMON/BLK51/XPTADJ(MAXVRT),XADJVT(2*MAXEDG).XEDGPB(MAXEDG) 
•COMMON/BLK52/  XLNKDW(2*MAXEDG).XLNKUP(2*MAXEDG),XLN1CED(2*MAXEDG) 
•COMMON/BLK53/  FOUND2 

•COMMON/BLK54/  XDEG(MAXVRT).XCARDE.XCARDV.XAVSAD 
•COMMON/BLK55/  EPROB(MAXEDG),RELB(MAXEDG),TOP.DISTGN 

•COMMON/BLK57/  TPTADJ(MAXVRT),TADJVT(2*MAXEDG).TEDGPB(MAXEDG),TM,TTOPB 
•COMMON/BLK58/  TLNKDW(2*MAXEDG).TLNKUP(2*MAXEDG),TLNKED(2*MAXEDG) 
•COMMON/BLK59/TDEG(MAXVRT),TCARDE,TCARDV,TCARDK,TAVSAD,TBRlDG(MAXEDG) 
•COMMON/BLK60/SEPVTX(MCARDV),TOPVTX 

•COMMON/BLK6 1  /  FORBEG(MAXEDG),BRlDGE(MAXEDG)JOPB,TOPFB 
•COMMON/BLK84/  VCPUJCPU 
•COMMON/BLK98/  DOMINT.NPIVOT 
•COMMON/BLK99/  COUNTE 

5.6.  Description  of  Subroutines 

Next,  the  subroutines  are  presented  and  briefly  described.  Subroutines  DELETE(V,PTR), 
SERIER(V),  and  DEG2R(V)  are  from  the  original  version  of  PolyChain. 

5.6.1.  SUBROUTINE  NONSEP(FNONSP) 

Description  This  subroutine  finds  the  nonseparable  components  when  the  graph  is  separ¬ 
able,  and  the  edges  which  removal  disconnects  the  graph. 

Input 


Output 


The  multilist  structure  and  the  logical  variable  FNONSP. 

A  list  containing  the  pointer  to  the  beginning  of  the  list  of  vertices  of  each 


component,  the  list  of  vertices  of  each  component,  and  the  list  of  bridges. 


5.6.2.  SUBROUTINE  XPUSH(VERTEX) 

Description  This  subroutine  puts  element  VERTEX  on  the  top  of  stack. 

Input  VERTEX  and  the  stack. 

Output  The  updated  stack. 

5.6.3.  SUBROUTINE  OUTSEP 

Description  Prints  the  output  listing  when  the  network  is  separable. 

Input  The  list  of  vertices  of  each  nonseparable  component. 

Output  The  vertices  of  each  nonseparable  component. 

5.6.4.  SUBROUTINE  FACTOR 

Description  This  subroutine  controls  the  basic  steps  of  the  factoring  algorithm. 

Input  The  multilist  structure  of  the  reduced  network  obtained  after  polygon-to- 

chain  reductions  were  performed. 

Output  The  K-terminal  network  reliability. 

5.6.5.  SUBROUTINE  REDUCE 

Description  This  subroutine  performs  series,  parallel,  and  degree-2  reductions. 

Input  The  multilist  structure. 

Output  The  updated  multilist  structure  after  ail  possible  simple  reductions  were  per¬ 


formed. 
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5.6.6.  SUBROUTINE  SER1ER(V) 

Description  This  subroutine  performs  a  series  reduction  on  vertex  V  not  in  set  K. 

Input  Vertex  V  and  the  multilist  structure. 

Output  The  updated  multilist  structure,  with  V  and  both  of  its  edges  deleted,  and 

with  a  new  edge  inserted.  This  new  edge  has  its  reliability  computed.  New 
cardinalities  of  V  and  E. 

5.6.7.  SUBROUTINE  DEG2R(V) 

Description  This  subroutine  performs  a  degree  2  reduction  on  vertex  V  in  set  K. 

Input  Vertex  V  and  the  multilist  structure. 

Output  The  updated  multilist  structure,  with  V  deleted,  along  with  both  of  its  edges, 

and  with  a  new  edge  inserted.  This  new  edge  has  its  reliability  computed. 
New  cardinalities  of  V  and  E.  The  updated  value  of  M. 

5.6.8.  SUBROUTINE  COPY 

Description  This  subroutine  copies  the  current  graph  for  later  use. 

Input  The  multilist  structure,  and  the  list  of  edges  that  are  bridges  for  the  current 

graph. 

Output  The  multilist  structure,  and  the  list  of  edges  that  are  bridges  for  the  current 

graph. 

5.6.9.  SUBROUTINE  SELECT(V1,V2) 

Description  This  subroutine  selects  an  edge  to  pivot. 

Input  The  multilist  structure,  and  the  edges  that  are  forbidden  to  be  chosen. 

The  nodes  that  form  the  edge  selected. 


Output 


5.6.10.  SUBROUTINE  CHKKGR(KGRAPH) 


Description  This  subroutine  checks  if  the  graph  is  a  K-graph. 

Input  The  multilist  structure,  and  the  list  of  vertices  of  each  nonseparable  com¬ 

ponent. 

Output  A  logical  variable  indicating  whether  the  graph  is  a  K-graph  or  not. 

5.6.11.  SUBROUTINE  GRAPHR(TOP) 

Description  This  subroutine  reconstructs  a  subgraph  of  the  computational  binary  tree  that 
is  pointed  to  by  TOP. 

Input  Pointer  TOP  and  the  multilist  structure  of  the  original  graph. 

Output  The  multilist  structure  of  the  reconstructed  subgraph. 

5.6.12.  SUBROUTINE  REMOVE(Vl,V2) 

Description  This  subroutine  removes  the  edge  incident  to  vertices  VI  and  V2  from  the 
subgraph. 

Input  Vertices  VI  and  V2.  The  multilist  structure. 

Output  The  updated  multilist  struaure,  with  the  desired  edge  removed. 

5.6.13.  SUBROUTINE  COLAPS(Vl,V2,) 

Description  This  subroutine  changes  the  subgraph  by  considering  the  probability  of  the 
edge  incident  to  vertices  VI  and  V2  as  being  equal  to  1. 

Input  Vertices  VI  and  V2.  The  multilist  structure. 


Output 


The  updated  multilist 


5.6.14.  SUBROUTINE  DELETE(V,PTR) 


Description  This  subroutine  deletes  the  element  pointed  to  by  PTR  from  vertex  V’s  adja¬ 
cent  vertices  list.  Three  cases  are  considered.  The  first,  when  the  element  is 
first  in  the  list.  The  second,  when  it  is  last  in  the  list.  The  last,  when  the  cle¬ 
ment  is  in  the  middle  of  the  list.  In  each  case,  the  element  is  deleted  by  a 
different  set  of  commands. 

Input  The  multilist  structure.  Vertex  V.  Pointer  PTR. 

Output  The  updated  multilist  structure  without  the  specified  element. 

5.6.15.  SUBROUTINE  FNDREL(REL) 

Description  This  subroutine  computes  the  reliability  of  a  reduced  subgraph. 

Input  The  stack  defined  by  EDGEVlf*),  EDGEV2,  and  DIRECT(*),  and  the 

current  reliability  REL. 

Output  Updated  r-sliability. 

5.6.16.  SUBROUTINE  FNDEDG(V1,V2,EDGE) 

Description  This  subroutine  finds  the  edge  incident  to  vertices  VI  and  V2. 

Input  Vertices  VI  and  V2.  The  multilist  structure. 

Output  The  edge  specified. 

5.6.17.  SUBROUTINE  OUTFAC(OUTPUT) 

Description  This  routine  prints  out  the  solution  after  the  factoring  algorithm  was  per¬ 
formed. 

Input  The  current  value  of  the  reliability  and  the  number  of  edges  selected  for 

pivoting. 


Output 


The  network's  reliability  and  the  number  of  edges  selected. 


6.  User  Manual 

The  user  manual  of  this  new  version  of  PolyChain  is  similar  to  the  manual  of  the  original 
version  [1].  In  this  section  we  first  present  a  guide  for  using  PolyChain  showing  the 
differences  when  using  the  VAX/UNIX  system  and  the  IBM/CMS  system.  The  input  file  and 
output  are  then  described,  and  a  test  problem  is  presented  to  illustrate  outputs  for  both  separ¬ 
able  and  nonseparable  cases. 

6.1.  Executing  Polychain 

Polychain  can  be  used  in  either  the  VAX/UNIX  system  or  IBM/CMS  system.  As  already 
mentioned,  only  I/O  related  code  is  system  dependent.  Therefore,  to  run  the  code,  first  the 
routine  that  gets  the  time,  date,  and  day  of  the  week  from  the  system  must  be  specified.  Then, 
the  dimension  parameters  and  the  COMMON  blocks  must  be  adjusted.  Finally,  an  input 
data  file  must  be  prepared.  These  three  topics  are  presented  below. 

6.1.1.  System  Routines 

The  first  step  in  running  Polychain  is  adjusting  the  code  to  run  in  the  desired  system,  either 
UNIX  or  CMS.  To  do  this,  the  suitable  system  routine  that  gets  the  time,  date,  and  day  of 
the  week  must  be  specified.  The  code  considers  both  possibilities,  so  that  is  just  a  question  of 
removing  or  adding  comments  to  the  lines  of  the  code  where  the  system  routines  appear, 
depending  on  which  one  you  need.  The  system  routines  are  described  below. 

For  IBM/CMS  use: 

MAIN  ROUTINE: 

CALLDATETM(DATTIM,23,VCPU,CTIME,TCPU) 

DATE  -  DATTIM(1:16) 

HOUR  -  DATTIM(  19:23) 

BEGINT  -  VCPU 


SUBROUTINE  OUTFAC,  OUTGRF  AND  OUTREL: 

CALLDATETM(DATT1M,23,VCPU,CTIME.TCPU) 
TMPCPU  =  VCPU  -  TMPCPU 
WR1TE(IOUT,300)  TMPCPU 


For  VAX/UNIX  use: 


MAIN  ROUTINE; 

CALL  FDATE(ERA) 
DATE  =  ERA(1:I0) 
HOUR  =  ERA(  12:20) 
YEAR  -  ERA(2l:30) 
CALL  DTIME(TIME) 


SUBROUTINE  OUTFAC,  OUTGRF  AND  OUTREL: 

CALL  DTIME(TIME) 
WRITE(IOUT,300)  TIME(l) 


6.1.2.  Dimension  Parameters 


The  second  step  in  running  Polychain  is  adjusting  the  dimension  parameters  MAXVRT  and 
MAXEDG  in  SUBROUTINE  INILST.  MAXVRT  is  the  maximum  number  of  vertices  and 
MAXEDG  is  the  maximum  number  of  edges  of  the  graph.  The  adjustment  of  these  variables 
is  needed  only  if  the  network's  dimensions  exceed  what  has  been  already  specified. 

After  adjusting  the  dimension  parameters,  all  COMMON  blocks  containing  arrays  must  be 
changed  accordingly.  Section  4.S  shows  how  the  arrays  must  be  changed. 


6.1.3.  Input  Files 


Inputing  data  in  PolyChain  is  very  simple  since  data  is  not  restricted  to  specific  columns  of 
the  input  line.  No  flag  is  needed  to  indicate  end-of-fiie.  The  first  line  of  the  input  file  con¬ 
tains  the  system  output  options.  One  value  must  be  entered  in  this  line  -  ECHOIN,  where. 


ECHOIN 


1  if  a  report  of  the  input  network  is  desired 
0  otherwise 


Next,  the  edges  are  specified,  one  in  each  line.  To  specify  an  edge,  enter  both  vertices  of  the 
edge  followed  by  the  edge’s  reliability.  The  numbering  of  the  vertices  should  be  sequential 


from  1  to  the  number  of  vertices  of  the  network.  If  a  vertex  is  a  K-vertex,  it  should  be  pre¬ 
ceded  by  a  minus  sign.  An  example  illustrating  an  input  file  is  given  below. 


Figure  5 


The  input  file  for  the  network  of  figure  5,  in  the  case  that  we  want  the  input  network  report  is 
given  below. 


-1  2  .5 

-1  3  .8 

2  3  .7 

2  -4  .6 

3  -4  .9 


6.1.4.  Program  Outputs 


In  this  section  a  test  problem  is  used  to  illustrate  the  program’s  output.  Consider  a  series- 
parallel  irreducible  network,  the  ARPA  computer  network,  in  figure  6.  The  reliabilities  (actu¬ 
ally  availabilities)  shown  in  figure  6  are  fictitious. 


PolyChain  generates  either  a  two  part  or  a  three  part  report  depending  on  whether  the  net¬ 
work  is  series-parallel  reducible  or  not,  respectively.  The  first  section  of  the  report  describes 
the  input  network,  edge  by  edge.  The  type  of  each  vertex  is  indicated,  K  for  K-vertex  and  nK 
for  non  K-vertex.  The  first  section  also  summarizes  the  input  network  data  and  core  usage. 
Network  density,  presented  in  the  first  section,  is  defined  to  be  the  ratio  of  the  number  of 
edges  of  the  input  network  to  the  number  of  edges  of  its  corresponding  complete  graph.  The 
second  section  of  the  report  indicated  whether  the  network  is  series-parallel  reducible  or 
irreducible.  This  section  contains  a  summary  of  the  reductions  performed  and  the  CPU  time 
before  the  beginning  of  the  factoring  algorithm.  In  case  the  network  is  series-parallel  irreduci¬ 
ble,  the  updated  value  of  A/  =  nn.  is  included  in  this  section  and  the  third  part  of  the  report 

j 

is  generated.  The  third  section  contains  the  K-terminal  network  reliability,  the  domination  of 
the  reduced  network,  the  number  of  pivots  performed,  and  the  CPU  time,  excluding  I/O.  The 
report  generated  by  PolyChain  for  the  above  file  follows. 


PCLYCHAIN  •  YERSICN  65.1 
PCLYCQN  TO  ChAIN  REOUCIICNS 
IK  KEThQRK  RELIABILITY 


INPIT  NEThCRK 
I  EDGE  I  1  VERTEX 


TYPE  11  VERTEX 


FACE  1 


DATE  :  FRI»  KCV  2S  1S65 
TIME  :  15:0E 


11 

19 

1 

NK 

II 

13 

1 

NK 

1  |.9aCCC0C0E'»00 

1 

•  1 

11 

13 

1 

KK 

11 

21 

1 

K 

I1.90CCC0Q0E400 

1 

1 

TYPE 


RELIABILITY 


.eQ€CCOCOE-*^OC 


•  BOCCCOOOE-^CO 


•BOCCCOCOE^OO 


.95CCC0COE-»00 


.SOCCCOOOE«00 


.90C0COC0E«C0 

41 

•BOCCCOCOE^OO 


.60CCC000E‘»0C 


•SCCOCOCOE^OO 


•9CCCC0C0E^0C 


•  9CCCC0C0E-»0C 


.9CCCC0C0E’»0C 


.90CCC0C0E«00 


•90CCCOCOE«00 


.90CCCOCOE«CO 


.90CCCOCOE-»O0 


.80CCC000E«0C 


.90CCC0C0E«0C 


•95CCC0C0E«0Q  1 


PCLVCHAIN  -  VERSION  85.1 
FCLVGON  10  Ct-Alh  REDtCIlCNS 
IK  NE1KQRK  RELIAEILITY 


INPtT  NETKERK 


PAGE  2 


DATE  t  FRI»  KCV  25  1585 
TINE  :  15:0E 


1  EDGE 

II 

VERTEX 

1 

TYPE 

II 

VERTEX 

1 

TYPE 

11  RELIABILITY 

1  22 

II 

9 

1 

KK 

II 

16 

1 

NK 

1  |.9CCCC0C0E«0C 

1 

1  23 

II 

16 

1 

NK 

11 

17 

1 

NK 

1 1.95CCC0C0E40C 

1 

1  25 

II 

■  A 

17 

1 

NK 

II 

»  A  A- 

lb 

1 

NK 

1  I.90CCCOCOE400 

1  25 

11 

■  A  A« 

18 

1 

NK 

II 

19 

1 

NK 

1  |.90C0COCOE«O0 

1 

26 

11 

15 

1 

NK 

II 

21 

1 

K 

1  |.95CCC0COE-»O0 

SIPNARY  OF  INPUT  NETWORK  DATA 


NIPEER  OF  VERTICES... 

NINCER  CF  EDGES . . 

NINBER  CF  K'VERTICES. 
NETWORK  DENSITY . . 


11 

16 

2 

0.125 


SUPNARY  OF  CORE  USAGE 


VARIABLE 

CURRENT 

USAC 

NAPE 

VALUE 

PAXEDG 

5000 

26 

PAXVRT 

2000 

11 

'•'vV 

» “j* 


•  *.*  V  .* 


.*/v 
' ,  *  .  •  •  • 


Kv.  V 
lA- 


T»V 


PCLYCHillN  -  VERSION  65.1 
PCLYCON  TO  CliAIN  REOOCTICNS 
IN  NETNORK  RELliBlLlTY 


DATE 

TINE 


FACE  3 


PRlt  KCV  2i  1965 
15:0E 


NEThORK  SERIES-PARALLEL  IRREDUCIBLE 
REDLCED  NETWORK 


I 


M  « 

1 

EDGE 

II 

VERTEX 

»  «»  • 

1 

TYPE 

II 

VERTEX 

1 

TYPE 

II  RELULILITY 

1 

1 

1 

IS 
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7.  Test  Problems 

Next,  the  results  obtained  by  PolyChain  applied  to  several  networks  are  given.  Some  of  the 
networks  are  obtained  through  a  random  network  generator,  while  other  tested  networks  are 
from  [5].  Problems  were  run  on  the  IBM  3081,  at  Berkeley.  The  code  was  compiled  on  the 
CMS  FORTVS  compiler  using  optimization  level  3.  CPU  times  were  measured  through  the 
DATETM  system  routine.  Table  I  contains  a  summary  of  the  networks  tested  and  table  II  a 
summary  of  lest  results.  Figure  7  shows  an  example  of  a  network  where  no  polygon-to-chain 
reduction  is  possible  for  any  set  K  chosen. 


Table  I  -  Test  Problems 
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Table  11  -  Test  Results 


Problem 

Edges 

%  Reduct 
Vertices 

on 

K-vertices 

Domination 

CPU  Time 

1 

69.2 

71.4 

0 

4 

0.12s 

2 

0 

0 

0 

6 

0.18s 

3 

0 

0 

0 

6 

0.18s 

4 

0 

0 

0 

6 

0.09s 

5 

0 

0 

0 

24 

0.72s 

6 

0 

0 

0 

24 

0.36s 

7 

0 

0 

0 

16 

0.48s 

8 

0 

0 

0 

11 

0.15s 

9 

0 

0 

0 

40 

1.40s 

10 

0 

0 

0 

448 

17.86s 

11 

0 

0 

0 

247 

2.38s 

12 

0 

0 

0 

11 

0.36s 

13 

0 

0 

0 

20 

0.17s 

14 

0 

0 

0 

11 

0.19s 

15 

50.9 

30 

25 

5063 

171.81s 

16 

43.3 

10 

0 

35 

1.10s 

17 

53.8 

46.7 

0 

52 

1.86s 

Figure  7  -  Ten  Vertex  Cubic  Graph 


8.  Conclusions  nnd  Recommendations 


This  report  discussed  the  design  and  implementation  of  two  features  that  enable  PolyChain  to 
treat  a  larger  class  of  problems.  The  implementation  of  both  features  maintain  the  charac¬ 
teristics  of  the  original  version  of  PolyChain  facilitating  further  extensions  and  enhancements. 


Further  testing  is  still  needed  to  ensure  the  code’s  correctness. 


To  insure  the  evaluation  of  the  K-terminal  network  reliability  in  a  more  efficient  form,  the 
program  should  apply  polygon-to-chain-reductions  in  addition  to  simple  reductions 
throughout  the  factoring  algorithm. 


In  the  case  of  separable  networks  a  code  using  PolyChain  as  a  subroutine  can  be  used  to  com¬ 
pute  the  reliability  of  each  nonseparable  component  and  then  compute  the  overall  reliability 
of  the  network. 
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